Initial spreading of low- viscosity drops on partially wetting surfaces 
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Liquid drops start spreading directly after brought into contact with a partial wetting substrate. 
Although this phenomenon involves a three-phase contact line, the spreading motion is very fast. 
We study the initial spreading dynamics of low- viscosity drops, using two complementary methods: 
Molecular Dynamics simulations and high-speed imaging. We access previously unexplored length- 
and time-scales, and provide a detailed picture on how the initial contact between the liquid drop 
and the solid is established. Both methods unambiguously point towards a spreading regime that 
is independent of wettability, with the contact radius growing as the square root of time. 
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How fast can a liquid drop spread over a surface? This 
basic question is relevant for applications ranging from 
printing and coating, to agricultural applications [J-Isl- 
In the final stage of drop spreading the dynamics are 
governed by Tanner's law, which relates the radius of the 
wetted area with time as r ~ £ 1/10 0,0- This extremely 
slow dynamics emerges from a balance between surface 
tension and viscous forces close to the contact line jjjj. 
Much less is known about the early stages of spreading, 
just after a spherical drop is brought into contact with 
a substrate at vanishing approach velocity. In contrast 
to Tanner's law, this dynamics is very fast cap- 
illary energy suddenly becomes available when the drop 
touches the solid, and this energy is concentrated into 
a singular point of contact. It has remained unclear 
whether or not the wetting conditions can influence such 
rapid inertial flows [13l4l5j. 

The initial stages of drop spreading are strongly rem- 
iniscent of the coalescence of two spherical drops of liq- 
uid, which very rapidly merge after contact is established 
16-21]. For low- viscosity liquids such as water, it is well- 
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known that the contact area between the drops grows 
as r ^ t 1 / 2 during coalescence. This can be explained 
from the balance of the inertial pressure inside the drop, 
~ p(dr I dt) 2 , and the capillary pressure, ~ jR/r 2 . Here 
p is the density, 7 the surface tension, and R the drop ra- 
dius. Interestingly, an identical scaling law was observed 
experimentally for water drops spreading on a completely 
wetting surface js|; apparently, the presence of a three- 
phase contact line does not affect the pressure balance 
during the initial phase of spreading. A rather different 
picture emerged, however, for drops spreading on par- 
tially wetting surfaces [9l4l l| . The dynamics was found 
to depend strongly on surface wettability, r ~ t a , with 
a non-universal exponent a that varies with the equi- 
librium contact angle [9|. This raises a number of in- 
triguing questions: How can the contact line, and the 
surface chemistry, affect the "coalescence" of a drop with 
a surface? Are the initial stages of spreading truly non- 
universal, or is there a hidden regime at smaller times? 
How is contact established on a molecular scale? 

In this Letter we reveal the initial spreading dynamics 
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FIG. 1: (Color online). Initial stages of drop spreading on 
partially wetting surfaces of varying wettability, (a) Molecu- 
lar Dynamics simulations of Lennard- Jones nanodrops (R = 
30nm), and (b) Experiments of water drops (R = 0.5mm). 
The top panels show side views of the liquid drop just after it 
has made contact with the partially wetting substrate. Lower 
panels are bottom views at times (a) t = 10, 35, 400 ps, and 
(b) t = 4,8,44 ps. The contact radius r(t) can be measured 
in time. 



of low- viscosity drops, using two complementary meth- 
ods: Molecular Dynamics simulations of Lennard- Jones 
nanodrops and high-speed imaging of experiments on 
millimeter-sized water drops [Fig. [T]. We access previ- 
ously unexplored length- and time-scales, and provide 
a detailed picture on how the initial contact between 
the liquid drop and the solid is established. While sim- 
ulations and experiments describe different dynamical 
regimes, both methods unambiguously point towards a 
universal spreading regime independent of wettability, 
consistent with the inertia-capillary balance r ~ t 1 ^ 2 . 
This contrasts the scenario proposed by [9] : At very early 
times after contact, the spreading exponent is indepen- 
dent of wettability, for contact angles ranging from com- 
plete wetting to very hydrophobic. 

Molecular Dynamics simulations. — The use of 
Molecular Dynamics simulations (MD) allows for study- 
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ing the initial contact between a liquid drop and a solid 
substrate down to molecular scale. To reveal the funda- 
mental mechanism of contact and subsequent spreading, 
we use a generic Lennard- Jones liquid. The advantage of 
the molecular approach is that, unlike continuum mod- 
eling, no assumptions on the moving contact line sin- 
gularity Q are needed. In MD, the wetting characteris- 
tics are directly controlled by the solid-liquid interaction, 
which determines the equilibrium contact angle 6 eq [22 1. 
The challenge, however, is to achieve sufficiently large 
drop sizes to recover a hydrodynamic regime. We there- 
fore study a quasi-two-dimensional geometry rather than 
axisymmetric drops, in which the system size in the y- 
direction is only 15 molecular sizes [cf. Fig. DIb)]. In- 
deed, contact problems such as coalescence are known to 
be essentially 2D phenomena and the same is expected 
here [H,[23|. 

We perform simulations on binary systems, in which 
two types of particles exist: fluid particles that can move 
around either in the gas or liquid phase, and solid parti- 
cles which are frozen on an fcc-lattice and constitute the 
solid substrate {24j|. All particle interactions are defined 
by the Lennard- Jones potential: 
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Here, is the interaction strength between particles i 
and j and aij the characteristic size of the atoms. This 
size is the same for all interactions, <r^ = a = 0.34 nm. 
The potential function is truncated at r c = 5<r (1.7 nm) 
where faj is practically zero. The mass of the atoms was 
set at 20 amu, and a timestep of 1.75 fs was found to be 
sufficient to accurately model these systems. The interac- 
tion strengths between the fluid atoms are en = 1.2fc#T, 
with ks the Boltzmann constant and T the temperature. 
The simulations are done in the NVT-ensemble, where 
the temperature is held at 300K using a thermostat, 
which is below the critical point for a Lennard- Jones fluid 
with the interaction strengths used. The fluid particles 
(amount: Ni = 304, 192) are initially positioned on an 
fcc-lattice (shaped with a cubic outline) far from the sub- 
strate (N s = 78, 300), but are free to move around and re- 
lax towards an equilibrium drop shape. Periodic bound- 
ary conditions are present in the lateral directions. The 
dimensions of the quasi-2D system are 240 nm, 5.1 nm 
and 120 nm in the x-, y- and z-directions respectively 
(Fig. [1]). The depth of the system is short enough to 
suppress the Ray leigh- Plateau instability and leads to an 
infinitely long cylindrical-cap shaped drop. 

The interaction strength between the solid and fluid 
defines the contact angle |22j |. We considered four differ- 
ent wettabilities (thus, four different values of e s i): e s i = 
(0.3, 0.4, 0.8, 1.2)k B T, giving eq = (115°, 100°, 60°, 0°) ± 
10°. The liquid density was measured to be p = 
664 kg/m 3 , and the surface tension 7 was measured in a 



separate, planar system at 7 = 0.017 J/m 2 . The viscos- 
ity of the liquid was measured in another separate system 
(Poiseuille-geometry) to be 77 = 3.64 • 10 -4 kg/(m-s). 
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FIG. 2: (Color online). Molecular Dynamics results, (a-d) 
"Bottom views" showing the molecules within 0.5nm from the 
substrate just prior to and after initial contact (t = —70 ps, 
—25 ps, —5 ps and 35 ps, respectively), (a) A small number 
of molecules from the vapor phase is close to the substrate. 
(b,c) Some fluctuating patches of higher density form, (d) 
A region of high liquid density nucleates at the substrate, 
from which one can measure r(t). (e) Radius of wetted area 
as a function of time for varying substrate wettabilities eq . 
Once the contact is established, we observe a power- law with 
exponent 1/2 for all values of eq . Results are displayed in Si- 
units on the right and top axis, and displayed in dimensionless 
form on the left and bottom axis. Contact radius r is rescaled 
with the initial dr op radiu s R, time is rescaled with the inertial 
time scale r p = \J pR 3 /j. 

The following procedure was used to bring the drop 
into contact with the substrate. First, the liquid is al- 
lowed to equilibrate far away (32 nm) from the substrate. 
During this stage, the drop will assume its cylindrical 
shape (R = 30 nm) and the liquid equilibrates with the 
vapour phase. Next, a body force is briefly applied on 
the fluid atoms until the drop moves towards the sub- 
strate. Just before the drop comes into contact with the 
substrate, the center-of-mass velocity of the drop is sub- 
tracted from the atom velocities such that the drop now 
'hovers' above the substrate. Due to the close proximity 
of the substrate (around 1 nm) the thermal fluctuations 
of the interface lead to first contact between the drop 
and the substrate after which the drop starts to spread. 
Using this method, the approach velocity of the drop to- 
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wards the substrate is zero and is not a parameter in this 
problem. 

What happens during the initial contact? In these 
early stages one cannot yet speak about a continuous 
liquid phase in contact with the solid. Instead, one first 
encounters the discrete, molecular nature of the fluid. 
Figures [2ja-d) show snapshots of the molecules that are 
within 0.5 nm from the substrate, represented as white 
dots. First a number of vapor molecules is randomly 
distributed over the surface [Fig. Ufa)]. As time pro- 
gresses, more molecules come into contact and form fluc- 
tuating "patches" of high liquid density at the substrate 
[Figs. Efb,c)]. The boundaries of these patches is ex- 
tracted by computing the number density field of atoms 
near the surface, and taking the iso-density contour half- 
way between the liquid and vapor density. Eventually, 
the patch becomes sufficiently large to span the entire 
depth of the quasi-2D simulation domain, from which we 
define the time of contact [Fig. 2(d)]. The exact defi- 
nition of t = does not influence our main conclusions 
below. From that moment, we track the boundaries of 
this wetting patch, which are the moving contact lines. 
The contact lines become sharper and are well-defined 
during the spreading, as can be seen in Fig. QJa). Note 
that while the average vapor density close to the surface 
is slightly larger than in the bulk, the surface coverage is 
very low and does not represent a precursor film. 

The key result of our MD simulations is that, once the 
liquid drop has established contact with the surface, the 
spreading follows a single power law. Figure [2fe) shows 
the contact radius r versus time on surfaces with varying 
wettability. For all contact angles eq we observe a scaling 
consistent with r ~ t 1 ! 2 (best-fit exponent: 0.48). It 
turns out that the contact angle of the substrate eq does 
have an influence on the spreading, but only through the 
prefactor: the exponent is always very close to 1/2. The 
prefactor increases as the contact angle decreases, such 
that drops spread faster on the more hydrophilic surfaces. 

Despite the very small length- and time-scales in these 
simulations, the spreading appears to be consistent with 
the hydrodynamic picture of inertia-dominated coales- 
cence The left axis and bottom axis in Fig. El^e) 
represent the data in dimensionless units, where lengths 
are scaled with the initial drop radius R and time with 
the inertial scale r p = yj pR? j^. In these units, the 
data span a range similar to previous experiments on 
millimeter-sized water drops [!, [s| . 

Experiments. — To verify whether the spreading be- 
havior observed in MD is also found experimentally, 
we carried out experiments in a previously unexplored 
regime. The required spatial and temporal resolution is 
achieved by high-speed recording of drop spreading from 
below, using transparent substrates and recording rates 
up to 600,000 frames/second. Typical images are shown 
in Fig. [TJd. The high-speed camera (Photron SA 1.1) is 
connected to a microscope (Zeiss Axiovert 25), which in 
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FIG. 3: (Color online). Experimental measurements of con- 
tact radius r plotted as function of time t for different sub- 
strate wettabilities. Results for three different equilibrium 
contact angles are plotted: clean glass (6 e ~ 0°; coated 
glass (Q e = 65°; A), teflon coated glass (9 e — 115°; 
The curves represent averaged data of repeated measurements 
(five or more for each eq ) per substrate for drops with radius 
R = 0.5 mm, showing the reproducibility of the experiments. 



combination with a 10X microscope objective (Zeiss A- 
plan, 10X) and reflective illumination gives a maximum 
resolving power of 2 /im/pixel. To capture a large pe- 
riod of the spreading process, frame rates used are in 
the range of 10-600kfps. Quere et al [8] and Bird et 
al [9], have shown that data for different drop sizes 
collapse by inertial rescaling. As here we focus on the 
influence of wettability, we consider only one drop radius 
R = 0.5 ±0.01 mm. 

To investigate the effect of wettability on the spread- 
ing, we performed experiments with water drops on three 
different substrates with different equilibrium contact an- 
gle eq : clean glass (almost perfectly wetting, eq « 0°), 
coated glass (0 r = 55°; a = 75°; eq = 65°) and teflon 
coated glass (0 r = 110°; a = 120°; eq = 115°). In or- 
der to avoid any condensation effects prior to spreading, 
the surrounding air is saturated with nitrogen gas. A 
thin needle is fixed at height D = 1 ± 0.02 mm above 
the substrate, thereby setting the initial radius R of the 
spreading drop (height D = 2R). With a syringe pump 
set at a constant volume rate of 1/il/min a pendant drop 
is grown at the needle tip, until it touches the substrate. 
This generates approach velocities < 2-10 -5 m/s, so that 
the outer gas dynamics has a negligible influence on the 
contact process. The radius of contact r(t) is determined 
from images as in Fig. [TJd, using a custom-made edge- 
detection algorithm in Matlab that finds the maximum 
image intensity slope in every frame. 

Our experiments confirm a single power law during the 
initial stages of contact. The measurements of the con- 
tact radius r(t) are shown in Fig. [3] on linear axes. One 




FIG. 4: (Color online). Experimental results. Contact radius 
r measured as a function of time (top and right axis) for three 
different equilibrium contact angles: clean glass (6 e « 0°; 
coated glass (Q e — 65°; A) 5 and teflon coated glass (6 e = 115°; 
©). The data shown, is an average of at least 5 measurements. 
The error bars denote the statistical error, which is larger than 
the measurement accuracy. On the left and bottom axis the 
data is normalized by the drop radius R and inertial time 
r p = yj pR 3 /~f respectively. A new regime is observed at ear- 
lier times, where the spreading is independent of the equilib- 
rium contact angle. The colored crosses are reprinted data by 
Bird et al. 0] (corresponding with: eq = 3°, x; 43°, X; and 
117°, x). The arrow indicates the smallest times that were 
accurately resolved in the study by Bird et al. 



observes that the data fall onto three different curves, 
corresponding to the three values of eq . The curves sep- 
arate about 0.1 ms after contact, showing a dependence 
on wettability at later times. However, the early-time dy- 
namics are independent of wettability. This is revealed 
in Fig. |U which shows the same data on log-log scale. 
We find that our data for different eq perfectly collapse 
at early times (t/r p < 0.04), and display an exponent 
close to 1/2 (best fit: 0.55). We reprinted the data by 
Bird et al. [9| (crosses) for completeness, and find a per- 
fect agreement with our data at t/r p > 0.1, which is the 
range of accurate resolution in ref. 9. The upper /right 
axis represent Si-units, while for the lower/left axis we 
employ the inertial scaling. Thus, the key point is that 
our measurements reveal a regime where wettability has 
no effect on the spreading at early times, not even in the 
prefactor. 

Discussion. — We have shown that early stage 
spreading of low-viscosity drops on a partially wetting 
substrate is independent of wettability. The wetted area 
is found to grow as r ^ t 1 / 2 , for all considered wetta- 
bilities: we find no influence on the spreading exponent 
by the presence of a contact line. This suggests that the 



mechanism of capillary wave generation, invoked to ex- 
plain # eg -dependent spreading exponents |9[, cannot be 
the dominant factor at very early times. Still, such cap- 
illary waves could be relevant for explaining the later 
stages of spreading in Fig. HJ where a departure from 
the 1/2 scaling is observed. However, this departure first 
arises when r/R > 0.2, in which case self-similarity of 
the bridge connecting the liquid drop and the substrate 
is lost, and scaling cannot be assumed a priori. It there- 
fore remains a challenge to explain the moment when the 
effect of eq becomes apparent in the experiments. There 
is, however, a subtle difference between the MD simula- 
tions and the experiments regarding the prefactor of the 
spreading law. For the experiments we observe a perfect 
collapse of the data on a single curve at early times. By 
contrast, the MD curves do not collapse, but the pref- 
actor increases with decreasing eq . We can point out 
at least two possible origins for this difference. First, the 
time- and length-scales of the two systems differ by orders 
of magnitude. In addition the simulations and experi- 
ments are not dynamically similar. While the rescaled 
results of Fig. [2] and |4] are very close, the Reynolds num- 
bers defined as Re = pr(dr / dt) /rj are very different: it is 
order unity in MD and order 100 in experiments. This 
suggests that the MD could be influenced by viscous ef- 
fects, and it would be interesting to further investigate 
spreading for highly viscous liquids 12, 2l[. Another key 
difference is the importance of thermal fluctuations at 
molecular scales. These are known to have a dramatic 
effect on the dynamics of drop pinch-off 25|, [26| , and it 
would be interesting to further explore their influence on 
spreading in the molecular simulations. 

Finally, although wettability does not affect the ini- 
tial rapid inertial flow in drop spreading, other cases are 
known to be strongly influenced by surface properties 



13l . Il5l . |27|. From a more general perspective, the com- 
bination of such inertial flows with a three-phase contact 
line therefore remains a challenge. 
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